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ABSTRACT 

There is evidence from radio-loud quasars to suggest that the distribution of 
the H/3 broad emission line (BEL) gas is arranged in a predominantly planar ori- 
entation, and this result may well also apply to radio-quiet quasars. This would 
imply that the observed full width at half maximum (FWHM) of the H/3 BELs 
is dependent on the orientation of the line of sight to the gas. If this view is 
correct then we propose that the FWHM can be used as a surrogate, in large 
samples, to determine the line of sight to the H/3 BELs in broad absorption line 
quasars (BALQSOs). The existence of broad UV absorption lines (BALs) means 
that the line of sight to BALQSOs must also pass through the BAL out-flowing 
gas. It is determined that there is a statistically significant excess of narrow line 
profiles in the SDSS DR7 archival spectra of low ionization broad absorption line 
quasars (LoBALQSOs), indicating that BAL gas flowing close to the equatorial 
plane does not commonly occur in these sources. We also find that the data is not 
well represented by random lines of sight to the BAL gas. Our best fit indicates 
two classes of LoBALQSOs, the majority (« 2/3) are polar outflows, that are 
responsible for the enhanced frequency of narrow line profiles, and the remainder 
are equatorial outflows. We further motivated the line of sight explanation of 
the narrow line excess in LoBALQSOs by considering the notion that the skewed 
distribution of line profiles is driven by an elevated Eddington ratio in BALQ- 
SOs. We constructed a variety of control samples comprised of nonLoBALQSOs 
matched to a de-reddened LoBALQSO sample in redshift, luminosity, black hole 
mass and Eddington ratio. It is demonstrated that the excess of narrow profiles 
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persists within the LoBALQSO sample relative to each of the control samples 
with no reduction of the statistical significance. Thus, we eliminate the possibil- 
ity that the excess narrow lines seen in the LoBALQSOs arise from an enhanced 
Eddington ratio. 

Subject headings: (galaxies:) quasars: absorption lines — galaxies: jets — (galax- 
ies:) quasars: general — accretion, accretion disks — black hole physics 



1. Introduction 



About 15% - 20% of quasars show broad UV absorption lines (originally defined as 
absorbing gas that is blue shifted at least 5,000 km/s relative to the QSO rest frame and 



displ a ying a spread in velo c ity of at least 2,000 km/si (IWevmann Ill997t iHewett and Foltz 



2003 



Reichard et al.l l2003t iTrump et all l2006t iGibson et al 1120091 ). Although evolutionary 



processes might be related to BAL outflows this does not seem to be the primary deter- 
minant. There is no indication that there is an excessive amount of radiation associated 
with reprocessed emission from dust as would be expected if the BALQSOs were in an evo- 
l utionary stage characterized by a large vo lume of dusty gas enshrouding the central AGN 



( IWillott et al.l 120031 : iGallagher et al 1120071). Also, the overall optical/UV spe ctra are strik- 



ingly similar to non-BALQSOs (jWeymann et al.l Il99ll : iReichard et al.l 120031 ) . Thus, it is 
widely believed that all or most radio quiet quasars have BAL flows, but the designation of 
a quasar as a BALQSO depends on whether the line of sight intersects the solid angle sub- 
tended by the outflow. This prevailing view is our fundamental assumption that motivates 
our search for a diagnostic that determines these preferred lines of sight. This is our primary 
working hypothesis and we explore possible effects that would be expected by preferred lines 
of sight. 

The "standard model of quasars" i s one of a hot ac cretion flow onto a black hole and 
a surrounding torus of molecular gas (lAntonuccil Il993l ). Theoretical treatments indicate 
that the BAL outflow can be an equatorial w ind driven from the outer regions of a luminous 
accretion disk that is viewed at low latitudes. iMurrav et all (119951) . or a bipolar flow launche d 
from the inner regions of the accretion flow ( Punsly 1999al b). Furthermore, Elvis ( 2000 ). 
proposed a purely phenomenological model in which the BAL wind begins near or at the 
accretion disk and flows at mid-range latitudes arranged so as to propagate just above 
the surface of the distant dusty torus. It has also been proposed that there is mo re than 
one w i nd source (polar and equa torial) for the BAL winds that c oexist in QSOs (iPunsly 
1999bl : iProga and Kallmanl 120041 ). It has been further argued by iBrotherton et all (120061 ) 
that observations rule out one preferred narrow range of lines of sight to the BAL wind. The 
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question is left open as to the distribution of lines of sights to BALQSOs, a question we hope 
to answer for LoBALQSOs in this study. 

There is very little direct evidence that we have on the line of sight to the B AL region. 
The most direct method implem ented so far is to use radio variability information (IZhou et al 
20061 ; iGhosh and Punslyl 120071 ). This information can be used to bound the size of the radio 
emitting gas then deduce that the radio emission must be viewed close to the polar axis 
and em anate from a relativis t ic jet, thereby avoiding the well known inverse Compton catas - 
trophe (IMarscher et all 1 19791 ; iKellermann fc Pauliny-Tothl Il969l ; iLind and Blandfordl 119851 ) . 
However, this method will only find sources with a polar orientation and it is limited to a 
subset of those sources that have sufficient radio flux to make the measurement and multiple 
epochs of observation with compatible sensitivity, a very small subsample. Thus, another 
orientation indicator is needed in this field of research. In radio loud quasars, the domi- 
nance of t he radio core relative to the large scale radio lobes is the standard orientation 
indicator (IWills and Browne Ill986l ). Core dominant objects are polar, lobe dominant ob- 
jects are viewe d at large angles from the jet axis. Almost all BALQSOs are unresolved 
with the VLA, iBecker et all (120001 ). but the radio fluxes are small and there is not sufficient 
sensitivity (and in many cases insufficient resolution) to get a meaningful estimate of core 
dominance. Furthermore, the spectra of BALQSOs are often similar to a special class of 
radio source, Gigah ertz peaked radio sources, for w hich this orientation indicator might not 
even be applicable ( Montenegro- Montes et all 120081 1 . 



A breakthrough paper in the field of BEL orientation was IWills and Browne I (119861 ) 
which showed that there was an anti-correlation between radio core dominance and the 
FWHM of H/3 in radio loud quasars that was explained by a broad line region (BLR) that 
was predominantly planar. The existence of a planar distribution of low ionization BEL 



gas se e ms to carry ove r to th e radio quiet regime as well. It was noted in lAntonucci et al 



(119891 ); iMaiolino et all ( 120021 ) that a planar distribution resolves the paradox that quasar 
Lya equivalent widths (and BEL equivalent widths in general) imply column densities that 
should show the Lyman continuum in absorption and with all profile s heavily damped. But , 
no quasar BEL cloud has ever been convincingly seen in absorption (lAntonucci et al 1119891 ). 
Since a large fraction of the solid angle around quasar must be covered with low-ionization 
BEL gas, yet no Lyman continuum absorption has ever been convincingly detected from 
the BEL clouds, it must be that quasars oriented favorably for BEL absorption are not 
seen as quasars. The only plausible way to do this is by placing the low ionization BEL 
clouds within the solid angle of the torus, so that these objects are not seen as quasars 
from our point of view. Furthermore, the large covering factor of the continuum by the 
equatorial planar BEL gas distribution creates an enormous volume of ionized BEL gas, 
sufficient to produce these large equivalent widths. Lines of sight far above the equatorial 
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plane (in the hole of the dusty torus, as expected for exposed quasar nuclei, lAntonuccil (119931 ) ) 
penetrate the relatively thin column densities orthogonal to the plane of the gas distribution, 
revealing the continuum in emission, with insufficient column density to produce damped 
Lyman absorption edges. This resolves the paradox. In conclusion, the notion of a planar 
distribution of low ioni zation BEL gas appears t o be applicable to all quasars and suggests 
that the main ideas of IWills and Browne I ( 119861 ) are relevant to predominantly radio quiet 
samples as well. There is no evidence of a difference between radio quiet and radio loud 
AGN in this regard. There are some differences in the spectra of radio loud and radio quiet 
quasars that have been n o ted in the literature that a re not entirely understood (see eg. 
Corbin and Francis I ( 119941 ); ICorbin and Boroson I ( 119961 )). but do not clearly have a bearing 
on this point 

Thusly motivated, we compiled the FWHM of low redshift, 0.4 < z < 0.8, LoBALQSOs 
in the SDSS DR7 archives to look for orientation information. This redshift range insures 
that both Mg II and H/3 emission lines are visible. High ionization, CIV, absorption will 
not appear in these spectra thus, our analysis is limited to LoBALQSOs. First we compare 
the quasars showing LoBALs with those that do not. Then we simulate the distribution of 
H/3 FWHM that would be expected for different lines of sight. This information is used to 
discuss the line of sight to the BLR of LoBALQSOs. 



2. The Distribution of H/3 in SDSS DR7 

In the DR7 release we found 102 LoBALQSOs (out of a sample of 10,069 quasars) in the 
redshift range, 0.4 < z < 0.8. The LoBALQSO determination is based on spectra with S/N 
> 7, using the criteria of continuou s absorption over a velocity interval greater than 1600 



km/s for a depth of at least 10% as in lZhang et all (120101 ) which see for a d etailed discu s sion o f 



this criteria. The fitting of the H/3 profile also follows the DR5 analysis of lZhang et all ( 120101 ). 
The distributions of H/3 FWHM of LoBALQSOs and non-LoBALQSOs from DR7 in this 
redshift range are compared in Figure 1. Because the bin populations of the LoBALQSOs 
are small, error bars were added to the bin frequencies. The horizontal error bars represent 
the width of the bins in Figure 1. The vertical (bin population) error bars are difficult to 
assess for small bin sizes, we chose to represent these with binomial statistics. It is notable 
from Figure 1 that there is a difference in the partial distributions of FWHM between the 
two samples for values smaller than 6,500 km/sec (a subsample that represents ~ 80% of 
the sources for both the total sample and the LoBALQSOs, separately). The LoBALQSOs 
are skewed toward lower FWHM in the partial distribution. The difference between these 
partial distributions is significant in both a K-S test (a null probability of 0.0178) and a 
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Fig. 1. — Comparison of the distributions of H/3 FWHM of LoBALQSOs and non- 
LoBALQSOs from DR7. Note the excess of narrow line profiles less than 2,500 km/s in 
the LoBALQSO population. The error bars are based on binomial statistics. 



-6 - 



Wilcoxon Rank Sum Test (a null probability of 0.0028). Thus, formally, we conclude that 
the distribution of LoBALQSO FWHM is poorly described by the distribution of FWHM of 
non-BALQSOs. To understand more than this, we consider the effects of a preferred line of 
sight to the BAL region and an elevated Eddington ratio in LoBALQSOs. 



3. Searching for an Alternative Physically Based Explanation 

Before adopting a line of sight interpretation of the excess of narrow emission lines in 
the LoBALQSO population in Figure 1, we explore the possibility that the effect is driven by 
a physical difference in the black hole accretion system between the LoBALQSOs and other 
quasars. We explore three potential physical parameters: the bolometric luminosity, Ly, 
the mass of the central black hole, and the Eddington ratio, T = Lboi/LEdd- To do this, 
we create "matched subsamples" of the SDSS DR7 nonLoBALQSOs that are matched in 
optical luminosity, black hole mass and Eddington ratio to LoBALQSOs. It is demonstrated 
that the difference in FWHM between LoBALQSOs and nonLoBALQSOs that was noted in 
Figure 1 is just as pronounced in the "matched" nonLoBALQSO samples. Obviously, all the 
quasar properties can not be identical between two samples in all moments of the probability 
distributions. One must choose the most important parameter and optimize for that one. 
We pay particularly close attention to subsamples that are matched with respect to T. 

In our opinion, L op ti ca i = ALa(5100A) is the most logical parameter to restrict if one is 
to construct matched subsamples for the following reasons: 

1. L op ti ca i is directly measurable and the other quantities are derived in consort with many 
assumptions. 



2. In IZhang et all (120101 ). Figure 6, it was shown that a comparison of 2MASS fluxes 
with SDSS fluxes that the rest frame near IR/optical colors are similarly distributed 
for LoBALQSOs and nonLoBALQSOs. Thus, this optically selected sample of LoB- 
ALQSO spectra show minimal reddening in the optical band, with a magnitude tha t 



is similar to that of the background galactic starlight contribution (IZhang et alll2010l ). 
We expect that L op u C ai is representative of the intrinsic optical luminosity up to a small 
correction on the order of the galactic starlight. With these small de-reddening cor- 
rections, Lopticai is an acceptable surrogate for L ho i in either the LoBALQSO sample or 
nonLoBALQSO sample 



Note that if one were to cull the sample of LoBALQSOs from an IR flux limited sample, we would not 
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3. It is shown in iGangulv et al I (2007) for HiBALQSOs and corroborated for LoBALQ- 
SOs in Figure 12 of Izhang et al ( 2010 ) that the biggest physical difference between 
BALQSOs and nonBALQSOs is that BALQSOs tend to have larger rest frame contin- 
uum fluxes. 



The significant difference in the rest frame continuum flux noted in point 3 above, is driven 
largely by the lack of low values of L optica [ within the BALQSO population. This is expected 
by dynamical considerations, low values of radiati on pressure do not produce enough force to 
drive outflows from the environs of a black hole ( Murray et all Il995l ; iPunslyl Il999bl ) . Figure 
12 of IZhang et all (120101 ) indicates that there are very few DR5 LoBALQSOs with L opt i ca i < 
4 x 10 44 ergs/sec. Thus, we construct matched samples by truncating the nonLoBALQSO 
DR7 sample on the low side. We create a variety of samples by varying the low end cutoff. 
We produced many such samples with a range of the derivable physical parameters designed 
to straddle those of the LoBALQSOs. Table 1 shows data for the various samples created 
with different cutoffs. 

In Table 1, we list the physical parameters of various subsamples for the sake for com- 
parison. Column 1 describes the types of sources, LoBALQSO or nonLoBALQSO, that 
comprise the sample. The second column gives the low end cutoff in L optica i that was used 
to define the sample. The next column is the number of sources in the sample. The fourth 
column is the resultant average L optical . The next two columns give the probability that the 
sources within the sample with FWHM < 6500 km/s are drawn from the same population 
as the LoBALQSOs with FWHM < 6500 km/s, for a Wilcoxon rank sum test and a K-S 
test. Columns 7 and 8 are the Eddington ratio and logarithm of the black hole mass, the 
derivation of which are given is the following paragraphs. 

The sample called, "LoBALQSO de- reddened" , is designed to compare the intrinsic 
luminosity of the LoBALQSOs to that of the nonLoBALQSOs. The reddening in the optical 
is small as noted above in point 2. Thus, our attempts to de- redden the sources are not 
likely to generate large errors, even if our assumptions are not highly accurate. Our main 
as sumption is based on the DR5 composite spectra for LoBALQSOs and nonLoBALQSOs 



m 



Zhang et all (120101 ). In i-band, it was found that the LoBALQSO composite continuum 



had about 0.25 to 0.33 magnitudes of attenuation relative to the nonLoBALQSO composite 
continuum. The wavelength of interest in the rest frame, associated with L opticah straddles 



expect this to be true. The SDSS optical selection criter ia will not find the more reddened LoBALQSOs 
that would be found in the IR samples. We also note that IZhang et all (|2010h showed that the LoBALQSO 
spectra are significantly reddened in the ultraviolet and therefore Ljjv = ALa(3000A) is not a good surrogate 
for L boi 
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Table 1: Physical Properties of Samples 



Sample 


cutoff 


Number 


Loptical 


P 


P 


r 


Log(M bh ) 


OIII EW 




10 44 c 




1Q 45c 


Wilcoxon 


K-S 




M Q 


A 


LBQSO a 





102 


1.05 ± 1.04 


1 


1 


0.274 ±0.247 


8.51 ±0.60 


19.0 ±29.4 


LBQSO a 





102 


1.32 ± 1.31 


1 


1 


0.307 ±0.276 


8.56 ±0.60 


19.0 ±29.4 


(de-reddened) 


















NLBQSO b 





10069 


0.67 ±0.75 


0.0028 


0.0178 


0.214 ±0.223 


8.44 ±0.46 


26.7 ±31.0 


NLBQSO b 


5 


4824 


1.03 ±0.95 


0.0008 


0.0054 


0.262 ±0.261 


8.57 ±0.44 


23.5 ±21.3 


NLBQSO b 


6 


3706 


1.18 ± 1.04 


0.0007 


0.0051 


0.277 ±0.271 


8.61 ±0.44 


23.2 ±21.0 


NLBQSO b 


7 


2919 


1.32 ± 1.13 


0.0005 


0.0049 


0.290 ±0.282 


8.64 ±0.44 


22.7 ±20.4 


NLBQSO b 


8 


2313 


1.47 ± 1.22 


0.0006 


0.0061 


0.304 ±0.293 


8.67 ±0.44 


22.3 ±20.1 


NLBQSO b 


9 


1860 


1.62 ± 1.32 


0.0005 


0.0055 


0.316 ±0.300 


8.70 ±0.43 


21.8 ± 19.6 


NLBQSO b 


10 


1478 


1.79 ± 1.43 


0.0005 


0.0078 


0.332 ±0.313 


8.72 ±0.43 


21.4± 19.3 


NLBQSO b 


9 d 


1551 


1.45 ±0.87 


0.0017 


0.0152 


0.321 ± 0.300 


8.66 ±0.42 


16.3 ±9.0 


NLBQSO b 


9 C 


1428 


1.41 ±0.68 


0.0029 


0.0258 


0.326 ±0.304 


8.64 ±0.42 


15.2 ±7.9 



a ow ionization broad absorption line quasars 
& Non-low ionization broad absorption line quasars 
c ergs/s 

d The O III line strength of the sample was cutoff from above to limit the sample to < 10 43 ergs/s 
e The O III line strength of the sample was cutoff from above to limit the sample to < 8 x 10 42 ergs/s 
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the high end of i-band and the low end of the less attenuated z-band. Thus, we expect 
that an average extinction of ~ 0.25 magnitudes of attenuation in the rest frame at 5100A 
is a good approximation. The observed flux is the sum of the reddened AGN spectrum 
and the contribution of the host galactic starlight. Therefore, on average, we acknowledge 
that the rest frame L opt u C ai computed from observed fluxes will be larger than it would have 
been if it were computed directly from the reddened continuum flux from the nucleus. The 
0.25 magnitudes of attenuation is therefore a liberal upper limit to the average value of the 
attenuation and the intrinsic AGN continuum flux will tend to be over-predicted from the 
rest frame L opt ucai computed from observed fluxes. As such, our estimates will tend to be 
an upper limit to the intrinsic average L opttica i of the LoBALQSOs. We used this number 
to de-redden L optica i: this amounts to multiplying each L op u ca i in the LoBALQSO sample by 
1.25 to get the intrinsic value of L optica i in the de-reddened sample. 

It is also of interest to compute the distributions of derived physical properties. First 
of all we estimate Mhh from the H/3 FWHM using the virial mass formula in equation (5) of 
Vestergaard and Peterson! (120061 ) . 



M bh (H(3) 

■^q6.91±0.02 



/ F(H0) V / AL A (5100A) \ 



V1000km/sy V. 10 44 ergs/s / 



0.50' 



(1) 



The logarithm of the values are tabulated in the eighth column of Table 1. 
reddened sample uses the de-reddened L optical in the formula. 



The de- 



The physical parameter of most reley ance in certain previo us studies is the Eddington 
ratio, T ( lBorosonll2002l ; ISulentic et alll2006l ). The interpretation of lBorosonl (120021 ) was largely 
driven by a principal component analysis of line properties. The BALQSOs tend to have 
weak, narrow H/3 emission lines, weak OIII lines and strong optical Fell lines, giving them 
low eigenvector 1 values. The claim is that a low value of eigenvector 1 is physical in origin 
and is driven by large a Eddington ratio (jBorosonl 120021 ). The corresponding explanation 
of the narrow FL3 excess that is seen in LoBALQSOs would be that these sources are very 
high Eddington ratio sources and this translates into a stronger ionizing flux which creates 
low ionization gas farther out in the gravitational potential (in units of gravitational radius), 
hence producing narrower H/3. 

The first step in the computation of T is to estimate L hn i from L optica i. Perhaps the 



most popular bolometric correction is that of iKaspi et all (120001 ) . L bot = 9AL A (5100). We 
note that this bolometric correction is predicated on L op u C ai representing the intrinsic optical 
luminosity. This is not the case for the LoBALQSOs. Hence, the need for the de-reddened 
sample. The resultant T values are tabulated in the seventh column of Table 1. Since this 
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parameter has received much attention in the past, in the context of LoBALQSOs, we look 
at a particular matched subsample of nonLoBALQSOs in detail. The sample with a cutoff 
at 9 x 10 44 ergs/sec is well matched to the de-redden distribution of T according to Table 1. 
We plot the distribution of T in Figure 2 for both the matched sample and the de-reddened 
sample. The distributions are indistinguishable in a statistical sense. The probability that 
the two samples are drawn from different distributions of T is only 0.249 according to a 
Wilcoxon rank sum test and 0.254 in a K-S test. 

Next, we plot the distribution of H/3 FWHM for the matched sample and the LoB- 
ALQSOs (reddening does not affect the LoBALQSO FWHM distribution) in Figure 3. The 
distribution is very similar to Figure 1. We conclude that T is a negligible physical factor 
with regards to producing the excess of narrow H/3 profiles for LoBALQSOs in the Figure 1 
and the LoBALQSOs in DR7 in general. 



In terms of the eigenvector space, iBorosonl ( 120021 ) claimed that BALQSOs are associated 



with very small eigenvector 2 (large accretion rate) and extremely small eigenvector 1 (large 
T). Table 1, shows that the matched subsample in Figures 2 and 3 has a slightly higher 
average accretion rate (using the intrinsic L optica i as a surrogate) and slightly higher T than 
the de-reddened LoBALQSOs. There are other samples in Table 1 that have slightly smaller 
values of L opt icai and/or T as well. For all these samples that surround the de- reddened 
LoBALQSOs in the 2-D parameter space, the excess of narrow lines in the de-reddened 
LoBALQSO sample persists with high statistical significance. Thus, there does not appear 
to be any physical properties associated with the location in the eigenvector parameter space 
that are responsible for the excess of narrow H/3 line widths in LoBALQSOs. 

Finally, we explore the relationship to eigenvector 1 in more detail. In the last col- 
umn of Table 1, we list the O III EW which has a large projection on eigenvector 1 



( jBoroson and Green! 119921 ) . Weak O III emission, in of itself, is not a plausible physical 
mechanism for launching a BALQSO wind, but it is possibly related to the BALQSO phe- 
nomenon on the basis of the previously quoted eigenvector space analysis (i.e., being very 
far from the accretion disk, the narrow line region is not a fundamental physical quantity 
that characterizes the black hole accretion system such as black hole mass, Eddington rate 
and bolometric luminosity). Exploring the luminosity of the O III emission can potentially 
provide insight into the interplay between eigenvector 1 and the existence of observed UV 
broad absorption lines. The LoBALQSOs in Table 1 have weak O III line strengths. None of 
the control samples that are created by means of a low end cutoff on L optica i have an average 
O III EW as small the LoBALQSO sample. We would like the control smaples to strad- 
dle the average value of the OIII EW in the LoBALQSO sample in order to elucidate the 
connection between the H/3 FWHM and the O III EW. Thus, we constructed two modified 
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Fig. 2. — Comparison of the distributions of T for LoBALQSOs and the "matched" sample 
of non-LoBALQSOs from DR7. The matched sub-sample is created from the full SDSS 
nonLoBALQSOS sample by the condition L optica i > 9 x 10 44 ergs/s. The error bars are 
based on binomial statistics. 
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Fig. 3. — The distributions of H/3 FWHM of LoBALQSOs and "matched" sample of non- 
LoBALQSOs from DR7. The matched sub-sample is created from the full SDSS nonLoB- 
ALQSOS sample by the condition L op u C ai > 9 x 10 44 ergs/s. Comparison to Figure 1 shows 
that the excess of narrow line profiles less than 2,500 km/s in the LoBALQSO population 
persists. The error bars are based on binomial statistics. 
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samples in the last two rows of Table 1, that are derived from the L op u C ai > 9 x 10 44 ergs/s 
sample with the additional constraint that the O III luminosity is bounded from above by 
Loin < 10 43 ergs/s or Lqiii < 8 x 10 42 ergs/s. A statistically significant correlation exists - 
although the strength of the correlation seems somewhat diminished relative to the sample 
(two rows above in Table 1) without a threshold on Loin- There is a trend in Table 1, 
as the Loin maximum threshold is lowered and the average Loin decreases in the three 
samples defined by L optica i > 9 x 10 44 ergs/s, the difference between the partial distribution 
of H/3 FWHM and the partial distribution of H/3 FWHM in the LoBALQSO sample becomes 
less statistically significant. This might b e expected on the basis of the correlation between 



H/3 FWHM and Loin in quasar spectra flBoroson and Green! Il992l ). However, the relevant 



conclusion drawn from these three samples in Table 1 is that the excess of narrow H/3 line 
profiles seen in LoBALQSOs is statistically significant more pronounced than would have 
been expected solely from the consideration of the correlation of Lqiii with H/3. 

In summary, we constructed samples of nonLoBALQSOs from SDSS DR7 that strad- 
dled the LoBALQSO SDSS DR7 sample in three physical properties, optical luminosity, 
central black hole mass and Eddington ratio. According to Table 1, regardless of whether 
the parameters were above or below those of the LoBALQSOs, there was always a highly 
statistically significant excess of narrow H/3 line profiles for LoBALQSOs. This was shown 
explicitly for a sample that was well matched in Eddington ratio. It is concluded that these 
parameters are not the physical origin of the excess narrow line profiles. This conclusion is 
not an artifact of reddening, as we verified this for de-redden LoBALQSOs. This further 
motivates our exploration of line of sight effects to produce the excess narrow lines seen in 
the DR7 LoBALQSO population. 



4. Simulating a Preferred Line of Sight 

In this section, we consider the distribution of non-LoBALQSO FWHM and simulate 
what the distribution would appear to be for different preferred ranges of the line of sight. 
This analysis will allow us to discuss the different orientation dependent models. Anticipating 
the use of the simulation to check "goodness of fit" of the binned data with a x 2 test, we 
want to minimize the number of free parameters in our model of the data. The purpose of 
this fit is to allow us to model the majority of the sources below 8,000 km/s which comprise 
92% of the entire sample. We consider the 8% of sources with larger FWHM as outliers 
that do not drive the orientation dependent properties of the other 92%. As such, we seek 
a simple single parameter fit that captures the skewness, peak and width of the distribution 
of FWHM below 8, 000 km/sec. A distribution with the desired skewness and exponential 
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type tail is the Gamma distribution that we parameterize as, 

/(FWHM) = Ar(FWHM)- (a - 1} [e- FWHM/y ] , (2) 

where /(FWHM) is the probability density, N is a normalization constant chosen to make 
the total probability equal to one and V is a free parameter that is in units of km/s. We 
eliminate the other free parameter, a, by choosing which fit minimizes the y 2 residuals. 
The result is indicated visually in Figure 4, the a = 5 plot is a more accurate fit than the 
a = 4 plot for the best fit values of V. This discriminant is driven by the large residuals 
generated at low FWHM for a = 4. Fine tuning the parameters does not change the 
fit significantly. No simple parametric distribution will fit the data in Figure 4 with high 
statistical probability in a \ 2 goodness of fit test. However, as stated above, this theoretical 
one parameter distribution with a = 5 and V a function of the distribution of the line of 
sight is adequate for estimating the gross properties of the skewness, variance and mean of 
the distribution that are induced by line of sight effects with the exception of the few percent 
of outliers residing at high FWHM. 

We note that the fit for the full sample of nonLoBALQSOs, GAMMA(q = 5, V = 775 
km/sec), works equally well for the matched subsample that was described in the last section. 
Figure 5 shows that the fit to the matched Eddington ratio subsample of nonLoBALQSOs 
from section 3 is very similar to that of the full sample of nonLoBALQSOs. The figure clearly 
shows very little difference between the bin populations of the matched and full samples. 
There is more scatter in the matched sample than for the full sample because there are less 
than 19% as many sources in the matched Eddington ratio sample. 

In order to simulate the change in the distribution of FWHM of DR7 QSOs in Figure 
4, if a preferred line of sight is specified, we consider the popular notion that the BEL gas 
is distributed in an equatorial pancake, orthogonal to the jet axis, with a random velocity, 
v r , superimposed on an equatorial velocity, v „, that is predominantly bulk motion from 



Keplerian rotation iJarvis and McClurd (120061 ); IWills and Browne I ( 119861 ) 



FWHM ~ 2Jv 2 r + vl sin 2 9. (3) 



In IWills and Browne I (119861 ) . it is suggested th at (v r /vp) 2 ~ 0-1. I n gene ral there must be 



a huge range of these ratios in QSOs, but the IWills and Browne I (119861 ) value represents 
a characteristic average value over the many QSO environments. Thus, we parameterize 
equation 3 as 

FWHM w VO.l + sin 2 9v b ir , (4) 

where vu r is a characteristic dispersion velocity of the BLR. The analytic forms in equations 
(2) and (4), allow us to transform the theoretical a = 5 distribution in Figure 4, to different 
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Distribution of SDSS Hp FWHM 

0.06 | 




HP FWHM (km/sec) 

♦ SDSS Data Gamma (5. 775) Gamma(4, 950) 



Fig. 4. — Comparison of the theoretical Gamma distributions of H/3 FWHM of non- 
LoBALQSOs from DR7 for GAMMA(o = 4, V = 950 km/sec) and GAMMA(a = 5, 
V = 775 km/sec). The distribution, GAMMA(a = 5, V = 775 km/sec), provides a much 
better fit to the sharply rising low side of the peak of the distribution. 
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Distribution of SDSS HP Matched SampleFWHM 

0.07 I 




Fig. 5. — Comparison of the theoretical Gamma distribution of H/3 FWHM of non- 
LoBALQSOs,GAMMA(a = 5, V = 775 km/sec) from Figure 4 and the distribution of 
H/3 FWHM of the matched subample of nonLoBALQSOs. The matched sub-sample is cre- 
ated from the full SDSS nonLoBALQSOS sample by the condition L optica i > 9 x 10 44 ergs/s. 
We also show the similarity of the distribution of the full sample of DR7 nonLoBALQSOs 
and the matched sample. 
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specified ranges of the line of sight. First, we i nvoke the r e sults of iBarthel I (119891 ) that 
are the basis of the "standard model of QSOs" in lAntonuccil (119931 ). QSOs are viewed at a 
random distribution of lines of sight within a range of about 45° (which represents the half 
opening angle of the molecular torus) of the jet axis (the normal to the accretion disk). The 
transformed distribution is determined by a change of coordinates defined by the mapping 
of the distribution of velocities, vu r to a specified line of sight by means of equation (4). 
This amounts to a change in the single parameter V in equation (2) (and the induced 
change in probability normalization, N). The distribution of vu r is defined by the parameter 
V = 1314km/s = Vblr- For a given distribution of line of sights, g(9), the distribution 
parameter V transforms as V — FVblr-, where 



/ g(9)V0.1 + sin 2 9 sin 9d 9 
jg{6) sin9d9 



(5) 



For the uniform distributions (in terms of solid angle) considered below, g{9) is just a step 
function. The results are presented in Figure 6. 

In Figure 6, the distributions have the following significance, 



The uniform distributi o n of l ine of sight angles, 0° < 9 < 45°, designated as "Uniform 
[0, 45]" is the IBarthel I (119891 ) empirical distribution for all QSOs 



The uniform distribut i on of l ine of sight angles, 0° < 9 < 25°, designated as "Uniform 
[0, 25]" is the IPunslvl Jl999bh theoretical distribution for LoBALQSOs 



The uniform di stribution of l i ne of sight angles, 80° < 9 < 85°, designated as "Uniform 
[80, 85]" is the iMurrav et all jl995h theoretical distribution for LoBALQSOs 



The uniform di stribu t ion of line of sight angles, 40° < 9 < 50°, designated as "Uniform 
[40, 50]" is the Elvis! (120001 ) phenomenological distribution for BALQSOs 



We discuss these simulated distributions in the context of the DR7 LoBALQSO FWHM 
distribution in the next section. 



5. Analysis of the Simulated Distributions 

In Table 2, we attempt to quantify how well the different theoretical models describe 
the data. The columns from left to right are first the model description, then the next 
two columns are the degrees of freedom and the x 2 statistic if the continuous distribution 
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Theoretical Distribution of Hp FWHM 



0.400 T 




18000 



♦ LoBALs Uniform [0, 45] Uniform [80, 85] Uniform [0, 25] 1 Uniform [40, 50] 



Fig. 6. — Comparison of the simulated distributions of H/3 FWHM for different distributions 
of the line of sight. Note that none of these represent the DR7 LoBALQSO data very well. 
As the line of sight becomes more polar the distribution is skewed towards smaller FWHM. 
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is binned in 1,000 km/s cells as in Figure I. In order to make a x 2 analysis plausible, 
we need a theoretically expected population of at least 5 per bin. Thus, we combine all the 
LoBALQSOs with a FWHM > 7,000 km/s into one final bin (a total of 7 bins are combined). 
The fourth column is the probability, P x 2, that we can reject the hypothesis that the data 
is described by the theoretical distribution. Note that we have abbreviated the "Uniform" 
Distribution by the symbol "U" in Table 2. We added the K-S "D" statistic in the fifth 
column which is a goodness of fit measure that is not sensitive to small bin populations. 
The K-S probability that we can reject the hypothesis that the data is described by the 
theoretical distribution is listed in the last column. 

The only fits that are not rejected by one or both of the statistical tests are the two 
component fits in rows 4 through 10 (see Figure 7). Any composite fit that is composed 
of the range 55% Uniform [0, 25] and 45% Uniform [80, 85] to 75% Uniform [0, 25] and 
25% Uniform [80, 85] can not be rejected on a statistical basis. Note that the equatorial 
model is rejected by this statistical analysis. The distribution, (0.67) (Uniform [0, 25]) + 
(0.33) (Uniform [80, 85]), represents the peak fairly accurately and it is also fits the high 
velocity tail reasonably well. 

6. Discussion 

There have been other discussions involving the possible orientation of the line of sight to 
the BAL gas which we discuss below. What distinguishes our methods from previous efforts 
is the very direct determination of the geometry. Previous discussions in the literature 
involved indirect physical arguments that invariably relied upon unverifiable assumptions. 
In this section we discuss our results in the context of previous discussions of orientation and 
line widths. 



Table 2: Goodness of Fit to the Line of Sight Models 



Model 


df 


x 2 




D 


Pks 


U [0, 45] 


6 


11.85 


0.934 


0.141 


0.97 


U [40, 50] 


6 


45.12 


> 0.9999 


0.235 


> 0.999 


U [80, 85] 


6 


226.79 


> 0.9999 


0.433 


> 0.999 


(0.75)(U [0, 25]) + (0.25)(U [80, 85]) 


4 


3.24 


0.481 


0.118 


0.87 


(0.70)(U [0, 25]) + (0.30)(U [80, 85]) 


4 


3.06 


0.452 


0.086 


0.71 


(0.67)(U [0, 25]) + (0.33)(U [80, 85]) 


4 


3.20 


0.473 


0.076 


0.66 


(0.60)(U [0, 25]) + (0.40)(U [80, 85]) 


4 


4.30 


0.633 


0.064 


0.60 


(0.55)(U [0, 25]) + (0.45)(U [80, 85]) 


4 


5.81 


0.876 


0.088 


0.73 


(0.55)(U [0, 25]) + (0.45)(U [40, 50]) 


4 


4.40 


0.645 


0.130 


0.93 
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LoBALQSO Distribution of Hp FWHM 
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♦ DR 7 LoBALs 



-0.67(Uniform [0, 25])+ 0.33(Uniform [80, 85]) 0.60(Uniform [0, 25])+ 0.40(Uniform [80, 85]) 



-0.55(Uniform [0, 25])+ 0.45(Uniform [80, 85]) 0.55(Uniform [0, 25]) +0.45(Uniform[40, 50]) 



Fig. 7. — Comparison of the simulated two component distributions of H/3 FWHM for 
lines of sight defined by, (0.67) (Uniform [0, 25]) + (0.33) (Uniform [80, 85]), (0.60) (Uniform 
[0, 25]) + (0.40) (Uniform [80, 85]), (0.55) (Uniform [0, 25]) + (0.45) (Uniform [80, 85]) and 
(0.55) (Uniform [0, 25]) + (0.45) (Uniform [40, 50]). The model (0.67) (Uniform [0, 25]) + 
(0.33) (Uniform [80, 85]) is tantamount to notion that 2/3 of the Lo B ALQS Os are polar 
wind sources that were proposed in the theoretical model of iPunsly Jl999bh and 1/3 of 
the LoBALQSos ar e equatorial winds sources proposed in the theoretical wind model of 
Murray et alJl995h . 
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6.1. Polarization 



Historically, the observed elevated polarization of BALQSOs has been used as an ar- 
gument that ordinary quasars app ear as BALQSOs because the line of s i ght is just above the 
dusty torus in the standard model ( Schmidt and Hines 1997 ; Ogle fl~997 ; Lamy and Hutsemekers 
20041 ). In axisymmetric geometries, equatorial lines of sight produce higher degrees of po- 
larization than polar lines of sight. Furthermore, the putative "equatorial" line of sight to 
the BAL gas is consistent with the idea that stray dusty gas above the torus is resp onsible 
for the reddening that is observed in the UV continuum (ISchmidt and Hines 1119971 ). This 
argument is indirect and very circumstantial. First of all, the elevated degree of polarization 
is not that high for most BALQSOs, the median polarization is > 1%. Theoretically, this 
value can be accommodated in any geometry due to modest attenuation of the con tinuum 
(see be low) by the BAL wind which enhances the prominence of the scatte red light (IPunslv 



999bJ) . Secondly, two o f the most p o 
(Il997h . and MRK 231, ISmith et al| ( 



are now 



arized BALQSOs, FIRST 1556+3517. iBrotherton et al 

13% in t 



1995) , have polarizatio ns 



know to be polar objects ( Ghosh and Punsly 20071; Revno 



more, both of these polar sources are heavily reddened (INajita et al 



dsetaj 



re nea r UV and 
2009h. Further - 



2000l : ISmith et~a1ll995h . 



Thus, these two well studied examples show that one can not use the enhanced polarization 
of BALQSOs and reddening of the continuum as an argument for equatorial BAL winds. 



Theoretically, it was pointed out in iPunslyl (jl999bl ) that the rare outlier polarizations of 
> 3% can be obtained in a polar line of sight if the assumption of axisymmetry is dropped. 
Modest reddening is a natural consequence of the polar wind model. There is significant 
attenuation from scattering in the t heoretical mod els as the continuum radiation field prop- 
he amount of attenuation is consistent 
30% to 50%. 



The UV must 



agates through the polar BAL wind (lPunslvlll9 99b) 
with the statistically based estimates of iGoodrich I (ll997f ) 
shine through larger optical depths (the inner regions of the wind) than the optical (outer 
regions of the polar wind), hence the polar models naturally produce reddening. Excessive 
reddening requires entrainment of dusty material on larger scales. 

The equatorial line of sight argument is very dependent on the assumption of axisymme- 
try for the scattering surface, which is apparently not justified (as noted above, the extreme 
outlier elevated polarizations can occur in polar geometries once the perfect axisymmetry 
assumption is dropped and evidenced by the high polarization in the polar sources MRK 
231 and FIRST 1556+3517). Within these equat orial models the Lo BALQSOs are viewed 
even closer to the equator than other BALQSOs ([Murray et all Il995h . However, the fits in 
Table 2 and Figure 6 indicates that this model cannot explain the DR7 data. The equatorial 
view enhances large FWHM and cannot explain the excess of small FWHM seen in the DR7 
LoBALQSO data. Figure 7 indicates that at most ~ 1/3 of the LoBALQSOs are equato- 
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rial outflows. The method used in this paper to determine the line of sight requires less 
physical interpretation of the data and less assumptions than the polarization and reddening 
argument. 



6.2. Radio Luminosity and Orientation 



The anti-cor relation b e tween radio luminosity and the degree of absorption in LoBALQ- 
SOs was used in iDai et all ( 120101 ) as an argument that most LoBALQSOs must be viewed 
through an equatorial BAL wind. They proposed that the LoBALQSOs tend to be weak 
radio sources because the radio emission is vastly weaker for perpendicular lines of sight to 
a relativistic jet. They made an anal ogy to powerful radio lou d sources in which radio cores 
are enhanced by Doppler beaming ( 1 Wills and Browne Ill986l ). It is far from obvious that 
there is a strong connection between radio loud quasar jets and the weak radio emission seen 
in most LoBALQSOs. The authors assume without justification that every LoBALQSO jet 
is identical (bulk Lorentz factor, power, etc.) and the difference in single point radio flux 
density is due to the line of sight. Core flu x is not a good orientation indicator, since the 
same type of logical argument presented in IDai et all ( 120101 ) leads to the erroneous conclu- 
sion that radio quiet quasars must be viewed near the equator because they are radio weak 
(compared to radio lou d sources). The tru e orien tation indicator is the ratio of core to lobe 
flux density as noted in IWills and Browne I ( 119861 ) . The data analysis presented in our paper 
suggests the alternative explanation that the mechanism that launches the BAL wind tends 
to inhibit the central engine that drives the powerful radio jet or the ability of that jet to 
propagate. 



6.3. Two Component Models 



The only oth er attempt that we know of to arr ive at the distribution of lines of sights to 
BALQSOs was in lBorguet and Hutsemekerd (J2010|). They fit CIV absorption (HiBALQSOs) 
line profiles in a parametric model. Unfortunately, the model is completely adhoc, it is 
comprised of a spherical fast wind with with a slower narrow equatorial wedge inserted. The 
bulk of the solid angle is called a polar wind, yet it is nearly spherical. Every parameter is 
adhoc, from the launch point, to the distribution of the broad emission line gas, the density 
and the vel ocity. They found the need for both components, in general, to reproduce the line 
shapes. In iProga and Kallmanl ( 120041 ). the numerical simulations showed that by adju sting 
the pa rameters one could get polar winds, equatorial winds, mid-latitude winds as in Elvis 
( l2000l ) or virtually any combination thereof. However, given the adhoc nature of the models 
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in 



Borguet and Hutsemekera (120101 ). it is very unclear if there is any physical significance to 



the parametric variations that led these authors to conclude that most sources are viewed 
along the equatorial plane. The driving piece of logic was that deep zero velocity absorption 
(P Cygni type) lines mean that one is viewing the BAL win d along the eqau t or. H owever, 
FIRST1556+3517 has deep zero velocity absorption in MgH. iBrotherton et all (Il997[ ). and it 
is polar. 



7. Conclusion 



In this paper, we used the distribution of H/3 FWHM in the SDSS DR7 data release 
to estimate the distribution of lines of sight to LoBALQSOs. Our analysis indicates that 
predominantly equatorial outflow is ruled out for LoBALQSOs. We also find that the data is 
not well represented by random lines of sight. The distribution has an excess of sources with 
narrow H/3 that is best fit by assuming two classes of LoBALQSOs, the majority (~ 2/3) 
are polar outflows and the remainder are equatorial outflows. By choosing a variety of 
subsamples of nonLoBALQSOs matched in redshift that straddle the values of luminosity, 
black hole mass and Eddington ratio of the de-reddened LoBALQSO sample, we find that 
the narrow line excess in the LoBALQSO sample persists in all cases and therefore eliminate 
the possibility that the excess narrow lines seen in LoBALQSOs arise from the physical 
properties of the central black hole accretion system such as Eddington ratio. 



We make one more comment on the important work of iGanguly et al I (120071 ). They 
measured the FWHM of the Mg II of HiBALQSOs and non-BALQSOs. There was no 
statistical difference. This implies random lines of sight to HiBALQSOs per the methods 
discussed here. However, they chose to combine mini-BALQSOs and BALQSOs to increase 
sample size. These types of sources are not BALQSOs in the conventiona l sense and the 
results might not be reliable. The DR5 statistical analysis of IZhang et all ( 120101 ) indicate 
that these types of sources (mini-BALQSOs have a large overla p in definition with the 
intermediate width absorption line sources of IZhang et all ( 120101 )) tend to resemble non- 
BALQSOs more than BALQSOs in many spectral properties. Furthermore, one must use 
caution in the interpretation of Mg II FWHM with HiBALQSOs, since there might be low 
level absorption in the continuum blueward of Mg II that might skew continuum estimates. 
It would be hard to segregate this effect and quantify its impact on the sample. This type of 
analysis is best done without resorting to resonance lines for FWHM estimates. Thus, H/3 
is preferred over Mg II. 



Based on these considerations mentioned above, it would be important to reproduce the 
analysis for LoBALQSOS presented here with a similar study for Hi BALQSOs. One could 
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select high redshift HiBALQSOs from SDSS and observe the H/3 profiles in the IR. 

We would like to thank Robert Antonucci for sharing his expertise and insights. We 
were also fortunate to benefit from many interesting comments from Paola Marziani. We 
would also like to thank an anonymous referee who offered many ideas that improved the 
manuscript. 
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